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Abstract. The intense interest in spin-based quantum information processing has 
caused an increasing overlap between two traditionally distinct disciplines, such as 
magnetic resonance and nanotechnology. In this work we discuss rigourous design 
guidelines to integrate microwave circuits with charge-sensitive nanostructures, and 
describe how to simulate such structures accurately and efficiently. We present a 
new design for an on-chip, broadband, nanoscale microwave line that optimizes the 
magnetic field driving a spin qubit, while minimizing the disturbance on a nearby 
charge sensor. This new structure was successfully employed in a single-spin qubit 
experiment, and shows that the simulations accurately predict the magnetic field values 
even at frequencies as high as 30 GHz. 



1. Introduction 

Magnetic resonance has been for over 60 years an essential tool to study the physical 
and chemical properties of matter pp. For most of its history, it has been a "bulk" 
technique [2J, where a large ensemble of spins is excited by means of radiofrequency 
of microwave excitations delivered by a centimetre-size resonant structure, designed 
according to microwave engineering guidelines |3j. The scope of application of magnetic 
resonance has evolved radically since the realization that both electron jl] and nuclear [5] 
spins in solid state can be used as qubits for quantum information processing [6]. For 
this purpose, the focus must shift towards the control and detection of individual spins 
at the nanometre scale. Such new and exciting line of research has driven the fusion of 
two traditionally separate disciplines such as magnetic resonance and nanotechnology. 

Applying high-frequency oscillating magnetic fields to a single spin is not any harder 
than to a large ensemble. The difficulty arises in integrating the microwave excitation 
with the ultra-sensitive detection techniques that must be employed to observe the signal 
of a single spin. Several examples of successful single-spin resonance in nanostructures 
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Figure 1: Optimized design for an on-chip broadband microwave transmission line for 
the control of a spin qubit. The inset shows a close-up of the nanoscale part of the 
structure, consisting of a short-circuit termination of the transmission line and a spin 
readout device. The dark sphere with arrow indicates the assumed location of the spin 
qubit. Alternative transmission line designs will be analyzed, while maintaining the 
same spin readout device structure. 

exist. They employed a wide variety of techniques to deliver the excitation to the spin, 
ranging from small coils [7J E] and waveguides [9] to nanofabricated coplanar striplines 
(CPS) [10] and coplanar waveguides (CPW) [UJ. These techniques tend to work well 
at relative low frequencies (~ 1 GHz) or when the spin under study and the method 
used to detect it are reasonably robust against heating and stray electric fields. A few 
groups [12, FTJ] have discussed the challenges they faced in integrating high-frequency 
(> 10 GHz) microwaves with single spins confined in GaAs quantum dots. 

The aim of our work is to explain and assess a set of useful design rules that can be 
applied to maximize the chance of success of high-frequency spin resonance experiments 
in nanostructures. We illustrate the pros and cons of different designs by taking as an 
example the architecture that was proposed for the control and readout of single donor 
spins [2], coupled to a silicon Single-Electron Transistor (SET) for spin readout [15]. We 
start with the broad microwave engineering framework to design a planar transmission 
line terminating in a loop close to the spin qubit. Then we present a novel on-chip 
transmission line for broadband localized B\ field generation (figure [I]). We describe 
a methodology for carrying out electro-magnetic field simulations on structures with 
large dimensional range, from millimetre down to nanometre size. We then use this 
tool to test the proposed transmission line and compare its performance with other 
designs. Finally, we present experimental measurements of our transmission line design 
and setup, and compare with simulations. 
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2. Design guidelines for the optimal planar loop 

The application we have in mind for the planar transmission line and loop is to drive the 
spin resonance of a single electron, while the spin state is read out through an energy- 
and spin-dependent tunnelling mechanism [TBI US]- The constraints posed by this type 
of application can be summarized as follows: 

1 - The electron spin Zeeman splitting Ez = h^B must be larger than the thermal 
broadening ~ 5k bT of an electron reservoir at temperature T. Here 7 = gfiB is the 
gyromagnetic ratio, g is the g-factor (g w 2 in Si, or g « —0.4 in GaAs), /xb is the Bohr 
Magneton, and B is an externally applied static magnetic field. The corresponding 
resonance frequency, i.e. the frequency that the transmission line must deliver to the 
nanostructure, becomes uq = Ez/h. For example, an electron temperature T ~ 200 mK 
would require a minimum resonance frequency v$ ~ 20 GHz. 

2 - The amplitude of the oscillating magnetic field B x produced by the loop should 
be maximized, to allow fast rotations of the spin. An electron spin will flip (i.e. undergo 
a 7r-rotation) in a time given by t n = h/(2gfisBi). The spin flip time t n should be much 
shorter than the dephasing time T£ to ensure high-fidelity spin rotations. is obtained 
from Ramsey or free induction decay experiments, yielding e.g. Trf ~ 50 ns in natural 
Si [T7]. Equivalently, Bi should be larger than the inhomogeneous width oesr of the 
electron spin resonance line, typically due to the coupling of the electron to the bath 
of surrounding nuclear spins. Typical values of ctesr are ~ 0.12 mT in natural Si [18J, 
and ~ 2.2 mT in GaAs [TU]. We note that, for the purpose of calculating the Bi value 
relevant to spin resonance, the rotating wave approximation usually holds. A linearly 
polarized oscillating field should be decomposed into two counter-rotating components, 
only one of which contributes to the spin rotation. Therefore, the B\ values obtained 
from a microwave simulator must be halved for the purpose of calculating e.g. a Rabi 
frequency. 

3 - Electric fields radiated from the transmission line should be minimized at the 
location of the spin qubit and readout device, since they can lead to unwanted effects 
such as photon-assisted tunnelling [10] , disrupt the operation of charge sensing devices 
(see section 4.3), and contribute to the local heating of the nanostructure. 

Given the constraints above, we will focus on simulations and design guidelines 
appropriate for spin resonance experiments in the 20 — 60 GHz range. However we 
note that a broadband planar loop that works well at v > 20 GHz will also exhibit 
good performance at MHz frequencies, since the losses and mode mismatches become 
less critical as the frequency is lowered. This is of relevance to the case where both 
electron spin resonance iy > 20 GHz) and nuclear magnetic resonance iy ~ 100 MHz) 



experiments can be performed on the same system, as in the case of dopant atoms [19 
2.1. Topologies of planar transmission lines 



We consider here two main topologies of planar transmission lines: coplanar stripline 
and coplanar waveguide. The fundamental difference between the two is that the CPS 
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Figure 2: Four planar loop designs, (a) A simple CPW on a PCB directly bonded to 
an on-chip CPS. (b) On-chip CPS with a CPW to CPS transition on PCB. (c) CPW 
matched to the termination, (d) Novel on-chip balun for maximized current at the loop. 



is a balanced line, meaning that it consists of two conductors, each having the same 
impedance to the surrounding ground planes. A CPS can carry microwaves in two 
modes: an odd mode, where the potentials of the conductors oscillate in opposite phase, 
and an even mode, where both conductors oscillate together with respect to the ground 
potential. The even mode is generally undesired, but can be excited in the presence of 
discontinuities and mismatches. The CPW instead is an unbalanced line, consisting of 
a single conductor, while the ground planes act as return lines. 

Coaxial cables - which we assume are going to be used to carry the microwave to 
the sample - are unbalanced transmission lines. Therefore we analyze the most typical 
situation in which a coaxial cable is first coupled to a CPW fabricated on a printed 
circuit board (PCB), then the microwave is further delivered to the chip containing the 
spin qubit. The CPS or CPW will have to be impedance- matched to the coaxial cable 
that delivers the microwave. The characteristic impedance of planar transmission lines 
are mainly a function of the width and thickness of the metal strips, the gap between 
coplanar strips and the thickness of the dielectric. Various transmission line analysis 
textbooks [201 EE] provide equations to calculate the impedance of many different planar 
transmission line types and configurations. 

Referring to the drawings in figure |2j we analyze four possible configurations of 
lines and loops. Where a CPS is employed, we describe how to realize the conversion 
from unbalanced to balanced mode, known as a balun. 



Figure 2a is the simplest solution, consisting of a transition between the CPW on 
the board to a CPS on the chip, realized by bond wires. This design can match the 
line impedances, but it does not provide a well-controlled unbalanced to balanced line 
transition. The CPS is terminated by a short circuit next to the spin qubit, with the 
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purpose of creating a node of the electric field - to minimize the disturbance to the 
SRD - and an anti-node of the current - to maximize the magnetic field that drives the 
qubit resonance. Here it's crucial to recognize that the desired node of the electric field 
only occurs for a perfect odd-mode transmission. Accidentally exciting an even mode 
will cause an anti-node of the electric field at the end of the line. This is a realistic 



danger for a CPW-CPS transition as crude that as shown in figure |2a| An on-chip CPS 
was used successfully by Koppens et al. [TU] to drive the magnetic resonance of a single 
spin in a GaAs quantum dot at v < 1 GHz, but contacting the CPS directly with a 
microwave probe instead of using bond wires. 

Figure 2b shows an improved design, where a CPW-CPS balun is first fabricated 
on the PCB, then the CPS on the board is bonded to the CPS on the chip. Here the 



conversion to a balanced line occurs in a more controlled way as compared to Fig. 2a 



and is less likely to excite an even mode of transmission. However, some risk still exists 
due to the use of bond wires to connect to the CPS on the chip. Planar baluns come in 
many forms and complexities. We chose the CPW to CPS transition presented by Chiou 
et al. [22] and analyzed by Mao et al. [23], which combines simple design, broadband 
operation and low insertion loss. 

Figure 2c shows a design based solely on CPWs, both on the PCB and on the chip. 
The CPW on chip terminates with short circuits on both sides. This type of design was 
used for instance by Fuchs et al. [IT] to drive the ESR of a nitrogen-vacancy centre in 
diamond at v — 0.49 GHz. Due to the absence of mode conversions, this design will 
have the widest range of operating frequencies and the lowest insertion loss. However, 
the nature of the CPW short implies that the current is divided equally amongst two 
loops, therefore only half of the signal received at the short can be exploited to generate 
a magnetic field at the qubit location. 

Finally, in figure 2d we propose a novel planar loop design, which consists of a 
short-circuited, on-chip CPW to CPS transition. In this design we keep the matched 
CPW structure at the interface between PCB and chip, and make the transition to CPS 
on the silicon chip. To achieve this, the transition needs to be scaled to accommodate 
the limited dimensions. 

In realistic experimental conditions, other parameters can influence the performance 
of the transmission lines. The chip will normally be mounted in a metallic enclosure for 
thermal and electromagnetic shielding. Delivering microwaves into the enclosure can 
excite cavity modes which can severely perturb the behaviour of the transmission lines. 
Therefore, the dimensions of the enclosure must be designed to minimize cavity modes. 
From microwave engineering textbooks [21], it is known that for CPW and CPS, keeping 
L/(S + 2W) > 1.75 and H/S > 2.5 ensures the enclosure will not affect the impedance 
of the transmission line by more than 1.5%; here L (H) is the width (height) of the 
enclosure, S and W are as defined in table [T} 

The shape and length of the wire bonds connecting the PCB and chip will also 
have an important effect on the signal losses at higher frequencies. Bond wires have an 
effective inductance which increases as they get longer and thinner. With the help of 
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Table 1: Dimensions used in all the planar transmission lines presented in this article. 
For CPW: S is the width of the centre metallic strip, W is the pitch between the signal 
strip and the ground plates. For CPS: S is the spacing between the balanced strips, W 
is the width of each metallic strip. 



Input Termination 
Type of TL S (fim) W (fim) S (fim) W (fim) 



CPW on PCB 


100 


60 


100 


60 


Balun on PCB 


100 


GO 


30 


300 


CPW on chip 


100 


GO 


1.6 


1 


CPS on chip 


30 


300 


1 


10 


Balun on chip 


100 


60 


7 


90 



simulations, we analyze this effect in more detail in section 3.3 



With regards to the design of the balun included in the topologies in figure 2b 



and figure 2d, we follow the method used by Mao et al. [23], where the CPW to CPS 
transition is divided into five sections. Each section tapers, crops and/or expands the 
transmission line to gradually make the transition keeping impedance variations small 
and fluid. The length of each section and the angle of the discontinuities needs to 
be optimized to minimize losses. When scaling our transitions to fit our PCB and 
chip dimensions, we ran a parametric sweep of simulations to optimize each transition 
section. The discontinuities at each section of the transition can excite the even mode 
of propagation, which can increase losses in the transmission line. This effect can be 
suppressed by connecting the ground planes with bond wires or air bridges where the 
discontinuities occur. 



2.2. Choice of simulation methods for tapered structures 

Modelling and simulation of the structures described in this paper is carried out using 
the Microwave Studio software package, from Computer Simulation Technology (CST- 
MWS). The choice of a suitable solver is based upon three requirements: 

1 - In most spin qubit device structures, the spin is buried at some depth below the 
surface on which the planar loop is fabricated. Therefore the electromagnetic simulation 
- the goal of which is to obtain the magnetic and electric field at the qubit location - 
must be carried out with a volume-based method. 

2 - We are interested in obtaining broadband results, since the operating frequency 
of the spin qubit might vary over a wide range, depending on the applied magnetic field. 

3 - All the structures described in figure [2] are sharply tapered, with edge lengths 
shrinking from millimetres to nanometres. The electromagnetic solver must handle 
correctly such wide range of structures sizes. 

The first requirement reduces our choice of available mesh types to hexahedral and 
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Figure 3: Tetrahedral mesh used for simulations. Mesh density is increased sequentially, 
adjusting to the sharp tapering of the structure. The rightmost inset shows the locations 
of the probes for all our simulation results: the magnetic field probe is at the donor site 
location (blue); the electric field probe is at the charge detector (red). 



tetrahedral. CST-MWS uses the finite integration method (FIM) [21] to support both 
these mesh types with time-domain and frequency- domain solvers. The chosen method 
should provide accurate results, in the least amount of simulation time. 

FIM in the time-domain is a leap frog algorithm similar to the finite-difference 
time-domain method [25]. It is useful for obtaining broadband results, as a single time 
domain simulation can be transformed to obtain an arbitrarily wide frequency response. 
However, this method requires a structured grid and only supports hexahedral meshing 
in CST-MWS. This means that for an accurate representation of our structure, the mesh 
around the small features will generate a great amount of unneeded mesh points at the 
large features. Simulation time for this scheme increases linearly on the number of mesh 
points and the minimum distance between mesh points. 

In the frequency-domain, FIM solves a set of Maxwell's equations for the entire 
volume, in one simulation. This makes the algorithm independent of the grid structure, 
which allows support for the more flexible tetrahedral meshing. Simulation times only 
depend on the total number of mesh points. The disadvantage of this scheme is that 
each simulation solves for only one frequency point, therefore it is necessary to run 
several simulations in order to obtain a broadband result. 

Attempting to accurately represent our tapered structure with hexahedrons will 
result in an unmanageably large mesh. Therefore the best option is to use the frequency 
domain solver with tetrahedral meshing. The immense difference in simulation times 
makes it worth running many simulations in the frequency domain, instead of one 
simulation with the time domain solver using hexahedral meshing. 

Once the solver and mesh type have been selected, the next key step is to 
produce a mesh that adequately represents all the materials, corners and edges of the 
structure, with the minimum number of mesh points. CST-MWS can be configured to 
automatically create an initial mesh and run a mesh adaptation algorithm. The software 
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will sequentially run a simulation and modify the mesh until a convergence is reached on 
scattering parameters (S-parameters) error. The algorithm adapts the mesh to increase 
its density at the electrically relevant locations (i.e. corners and edges). However, the 
automated meshing scheme fails to accurately adapt the mesh to accommodate the 
large feature size range in our sharply tapered structure. It is therefore necessary to 
manually divide the model into different sections and assign different meshing densities 
to each section (see figure [3]). With the mesh adjusted to both feature size and electrical 
relevance, we can obtain accurate results with minimal mesh points. 

Boundary conditions must be set for each of the planes that define the limits of 
the modelling space. Boundary conditions can be defined to either minimize reflections 
(open space), or behave like a perfect metallic object (electric). To simplify simulations, 
we model a smaller PCB than what would typically be used in an experimental setup, 
and we set the boundary conditions to open space on every side of the structure. This 



gives accurate results as long the guidelines for enclosure size discussed in section [2TT 
are followed. 

Electric and magnetic field probes can be placed anywhere in the three dimensional 
space of the structure. As shown in figure |3j the electric field probe is placed directly 
underneath the gate of the charge detector; while the magnetic field probe is placed 
at the donor site (i.e. where the spin to be controlled is located). We extract only 
the component of oscillating magnetic field that is perpendicular to the static magnetic 
field. 

All the models presented in this article comprise a modelled coaxial input port, 
a PCB with a planar transmission line, and a Si chip with a planar loop next to a 
spin readout device (SRD) similar to the one used by Morello et al. [TH]. We always 
assume 1 mW (0 dBm) of power at the coaxial input port of the simulated structure. 
We modelled the PCB using the characteristics of the Rogers RO3010 laminate with a 
thickness of 640 //m and 35 /iin copper (Cu) cladding. The planar transmission line on 
the PCB bridges the coaxial port and the chip. It is connected via bond wires to the 
planar loop, which is assumed to be fabricated in 100 nm thick aluminium (Al), on top 
of a silicon (Si) chip with a surface of 1.2x1.2 mm 2 and a thickness of 500 /im. As shown 
in Figure [TJ the loop is located at the centre of the chip surface, with a minimum width 
of 100 nm at the short circuit. The SRD is placed approximately 130 nm away from the 
loop. The conductivity of the Cu and Al lines is increased by two orders of magnitude 
as compared to the room-temperature textbook values, to account for the use of the 
device at cryogenic temperatures [26J [27J [28] ■ We assume the structure is always placed 
in a magnetic field large enough to suppress the superconductivity in Al. The bonding 
pads of all Al strips have a minimum area of 100 x 100 firn 2 . Unless otherwise stated, the 
bond wires have 25 /xm diameter and a length of approximately 200 /im. Table [T] shows 
the relevant dimensions of all the planar transmission line types we present. Weiland et 
al. [29] provides a more detailed explanation of the simulation parameters described in 
this section. 
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Figure 4: Characterization of the planar transmission lines, (a) 2D plot of the 
perpendicular component of the magnetic field amplitude perpendicular to the surface of 
the chip, (b) 2D plot of the absolute value of the electric field, (c) Magnetic field decay 
with distance from the loop. Lines in this plot compare simulations of the structure 

with (— — ) and without ( ) the SRD on the chip. The approximate location of the 

SRD and the spin qubit site are shown. The fields in all figures are calculated at 50 
GHz. 



3. Simulation Results 

3.1. Electric and magnetic fields along the transmission line 

Using the simulation setup described in the previous section, we started by carrying out 
a characterization of the electric and magnetic field profiles along the short circuited 
coplanar transmission line shown in Figure [Tj Figure 4a shows the amplitude of the 



perpendicular component of the magnetic field Bi generated around the loop. B x clearly 
decays with distance from the short, which shows the importance of placing the loop as 



close as possible to the spin qubit. In figure 4c we plot the decay of Bi as a function of 



distance from the loop, comparing the case with and without the spin readout device 
next to the loop. The SRD has a screening effect that further decreases the value of B\ 
beyond its location. 

The voltage difference between signal and ground lines generates an electric field, 



which decreases to zero (figure 4b) as it approaches the voltage node at the short-circuit. 



However, the resonant stub behaviour of the on-chip line (see section 3.2) causes an 
additional electric field to be emitted at the end of the line. This electric field can be 
responsible for a degradation of the performance of the SRD, and heating of the electron 
layer nearby. 

The surface current monitor allows us to estimate the current flowing through the 
loop. We find that 200 fiA of current flowing through the planar loop will generate 
0.1 niT of magnetic field at the spin location. 

The skin depth effect can reduce the effective conductivity of transmission lines. 
The skin depth 5 can be calculated as 5 = \j2p/ujp where p is the resistivity of the 
conductor, uj is the angular frequency of the current and /i is the absolute permeability 
of the conductor. Between 30 and 60 GHz, the skin depth of an aluminium line at 
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cryogenic temperatures (p = 269 pfi/m) is 33.5 nm to 47.2 nm. The smallest cross- 
section of the loop is 100x100 nm 2 , which therefore still allows for maximum current 
flow at the short. 



3.2. Performance comparison between CPS/CPW topologies 



The four planar transmission line models presented in section 2.1 were simulated and 
the results are compared in figure [5j The Sn parameter, which represents the reflected 
power, is a useful measure of the broadband performance of the lines. An ideal 
short-circuited transmission line should be fully reflective, yielding Sn = dB at all 
frequencies. Sn < dB indicates radiative losses, and we indeed observed that the 
behaviour of Sn correlates with the electric field radiated by the line towards the SRD. 

The simulations show that the fully-CPW structure shown in figure 2c is the closest 
to having the ideal behaviour Sn ~ dB from to 60 GHz and, accordingly, the smallest 
electric field radiated at the SRD. This is not surprising, since this design is matched in 
impedance and propagation throughout. However, as mentioned before, this comes at 
the price of roughly halving the strength of B\ available at the spin qubit location. 

The on-chip balun design of figure 2d generates the highest B\ fields through most 
of the spectrum. It radiates adequately low electric fields, comparable to the CPW 
at low frequencies, but shows a resonance at around 59 GHz, where the electric field 
increases accordingly. This phenomenon, known as "stub resonance" , is caused by the 
discontinuity between the PCB and the Si chip, which makes the on-chip section of the 
transmission line behave like a short-circuited resonator [3J. A short-circuited line acts 
as a parallel R-L-C circuit when its length equals a quarter wavelength (A/4). In this 
type of resonator, a A/4 standing wave is formed by having a minimum of the current 
at the input and a maximum at the loop. The on-chip line modelled here has a length 
of 550 pm, thus the resonance should be observed at A = 2.2 mm. Frequency and 
wavelength in microwaves relate to each other through the propagation velocity, given 
by v — d I y/^ejfi where e e ff is the effective dielectric constant of the medium through 
which the wave is travelling. Due to the fact that waves in planar lines travel at the 
interface of two different media (in this case vacuum, e r = 1, and silicon, e r = 11.9), 
obtaining e e // for our on-chip balun is non-trivial. Assuming an ideal planar line, half 
of its field travels in the substrate (filling factor of 50%), and its effective dielectric 
constant is the average between the two media. With this assumption we can estimate 
e e ff ~ 6.45, and a stub resonance at / = v/X — 53.7 GHz. The stub resonance observed 
in the simulations is shifted to a higher frequency due to the low inductance of the bond 
wires interfacing the PCB and chip. In section 3J3 these effects are analyzed in further 
detail. 



The design with propagation and mode mismatch (figure 2a) shows the poorest 
performance at high frequencies, radiating more than two times the amount of electric 
field compared to other designs. 





Figure 5: Simulation results comparing the frequency response of four planar 
loop designs: On-chip CPS with no transitions ( — • — ); On-chip CPS with PCB 
balun ( ); On-chip CPW ( ); On-chip balun ( ). (a) Sn parameter. 

(b) Magnitude of the perpendicular component of the magnetic field at the donor site. 

(c) Absolute electric field at the charge detector, in this case the SET island. 



Nanoscale broadband transmission lines for spin qubit control 



12 




Figure 6: (a) Comparison of Sn parameters for different bond wire lengths and types: 

200 /im bond wire (— — ); 500 /xm bond wire ( ); 150 fim bond wire ( ); 

200 fim ribbon bond 70 fim wide by 10 /im thick ( — • — ). (b) Electric and magnetic 

fields obtained from the model with (— — ) and without ( ) the bond wires bridging 

the ground plates of the on-chip CPW. (c) Comparison between: the standard model 
with the qubit located outside the loop (— — ); model with the transmission line elevated 

with calixarene, with the qubit located inside the loop ( ); standard model with 

a thin film of calixarene surrounding the transmission line and gates ( — • — ). The 



coplanar structure under study on all subfigures is the on-chip balun shown in figure 2d 



The device and probe locations with respect to the substrate are identical in all models. 



3.3. Further optimization 



As explained in section 32, the discontinuity between the PCB and the chip causes a 
stub resonance at high frequency Bond wires have an effective inductance that increases 
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with length and decreases with thickness |30| . In figure [6^, we show that with very long 



bond wires the resonance is very pronounced at the frequency calculated in section 3.2 
Shorter bond wires push the resonance to higher frequencies. This behaviour suggests 
that the bond wires can be modelled as an inductance in parallel with the R-L-C 
equivalent circuit of the on-chip resonating stub. Therefore, minimizing the effective 
inductance from the bond wires helps removing the resonant effect. 

Using very short bond wires can completely remove the resonant effect, but bond 
lengths of less than 200 /zm can be challenging to realize. A more viable solution is to use 
ribbon bonds, which have a much smaller effective inductance than standard bond wires. 



We have seen in figure [4b] that the electric field is mainly radiated towards the 
outside of the loop. Therefore, placing the spin readout device on the inside of the loop 
should decrease the electric field observed by the qubit, while maintaining high oscillat- 
ing magnetic fields. In the experiment performed by Koppens et al. [10], the SRD was 
indeed positioned on the inside of the loop, by separating the transmission line and the 
device with a thin film of dielectric. We can adjust our model to show results for this 
setup. We model a 100 nm coating of a material with e e ff = 7.1 (i.e. calixarene [31J) on 
top of the Si substrate and lay the planar loop on top, with the magnetic probe 50 nm 
inside the loop. We show in figure [6]d the simulated B 1 and electric field amplitudes 
of the shorted on-chip balun with and without the dielectric. The results suggest that 
using the dielectric layer to place the device inside the loop does improve the magnetic 
to electric field ratio at frequencies below 1 GHz and above 40 GHz. We also find that 
by adding the dielectric alone, without modifying the position of the loop, a shielding 
effect causes a reduction of the electric field at the SRD, without affecting the magnetic 
field. 



We explained at the end of section 2.1 the importance of bridging the ground 



plates of a CPW at discontinuities, in order to suppress the excitation of even modes of 
propagation. The simulation results in figure ^jp show how removing the bond wires that 
bridge the ground planes lead to a significant deterioration of the frequency response. 

4. Benchmark with experiments 

We have fabricated and tested a complete spin qubit structure as sketched in figure [2] and 



shown in figure 7a Due to the difficulty in making wire bonds of such short lengths, we 



were not able to bridge the ground plates of the on-chip balun in our fabricated model. 
4-1- Stub resonance of the on-chip CPW /CPS balun 

As a first experimental benchmark of the reliability of the models discussed above, we 
made an Sn measurement using a 10 MHz to 50 GHz Agilent PNA microwave network 
analyzer. Instead of bonding the planar line to a PCB, we have used a microwave probe 
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Figure 7: (a) Scanning electron microscope image of a spin readout device coupled to a 



microwave transmission line with on-chip balun as sketched in figure 2d Coloured dots 
show microwave probe locations corresponding to each Sn result, (b) Simulated (solid 
lines) and experimental (dots) Sn parameters for the on-chip balun. 
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Figure 8: SET current as a function of gate voltage, in the presence of a continuous- wave 
microwave excitation at 30 GHz. The data was taken on the same device used by Pla 
et al. [19]. Notice the broadening of the Coulomb peaks at high microwave power The 
microwave powers quoted in the figure are the output powers at the source. The power 
at the chip is approximately 30 dB lower due to loss along the coaxial line. The electric 
field radiated by the loop is extracted from the width of each of the peaks labelled on 
the figure. 



station, which terminates in a 3-terminal probe that can be placed directly on the chip 
at various locations. This allows us to study the stub resonance of the on-chip balun 
structure, by varying the distance between the edge of the chip and the contact point 
of the probe. To simulate this setup accurately, we removed the PCB structure from 
the model shown in figure 2d and connected bond wires from the coaxial port to the 
on-chip line, at three different microwave probe locations. Figure 7b shows a very good 
agreement between the measured and simulated Sn values. 
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4-2. Oscillating magnetic field at the spin qubit 

The on-chip balun planar loop design shown in figure [2] has been successfully employed 
to demonstrate coherent operation of a single-donor spin qubit in silicon [19]. The 
power-dependence of the Rabi oscillation frequency yielded an oscillating magnetic field 
at the donor B\ = 43 /xT at 30 GHz and dBm power at the source. 

In these experiments, the loss of the coaxial line connecting the source to the 
device was ~ 30 dB at 30 GHz, measured independently. Rescaling the simulations 
results in figure to a power of -30 dBm at the PCB and accounting for the rotating 
wave approximation (section |2|), our simulations predict 50 /iT of B\ at the donor 



(figure 5b). We consider this to be a very good agreement between the simulations and 
the experiment, and highlight the fact that the microwave simulation is quite successful 
in predicting surface currents and magnetic fields. 

4-3. Electric field radiated to the spin readout device 

The use of an electrostatically-induced SET [32] as the spin readout device allows us to 
estimate the magnitude of the electric field produced by the planar loop at the location of 
the SET island [33J . The current through the SET as a function of gate voltage exhibits 
characteristic Coulomb peaks. The width of these peaks in units of gate voltage AV g is 
a function of temperature, source-drain voltage bias and transparency of the potential 
barriers around the SET island. Because the SET island has a floating potential, an 
oscillating electric field Ei produced by the loop will contribute an additional Coulomb 
peak broadening AV 9 e- A simple parallel-plate capacitor model between the gate and 
the SET island yields the relation Ei = AV g E/d, where d is the distance between gate 
and SET island. In this case, d ~ 8 nm represents the thickness of the insulating Si02 
layer. 

A scan of the SET current as a function of the gate voltage shows broadening of the 



by writing AV gE = yj AV* T - AV g % 
peak with (without) the microwave 



Coulomb peaks when a microwave excitation is applied to the loop (figure 8a). The peak 
broadening AV 9 e can be extracted from figure 8a 
where AV g x (AV^o) is the width of the Coulomb" 



excitation. This calculation yields E\ « 380 kV/m from the data in figure 8a, taken at 
30 GHz and dBm power at the source. 

Rescaling the results in figure [6b for a microwave source to device loss of ~ 30 dB 



(see section 4.2), we obtain a prediction of E\ ~ 20 kV/m. This large discrepancy 
suggests that there could be additional mechanisms causing Coulomb peak broadening 
when the microwave excitation is applied, and/or that the simulation does not capture 
all the factors that influence the electric field patterns. For instance, additional thermal 
broadening of the Coulomb peaks could arise from radiated heat due to the sizeable 
current flowing through the loop. Additional electric fields could arise from cavity 
modes excited in the metallic enclosure used in this experimental setup. 
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5. Conclusions 

In conclusion, we have presented an extensive study on the integration of microwave 
transmission lines in nanostructures, for the purpose of optimizing the design of magnetic 
resonance experiments aimed at single spins. 

We have discussed a variety of topologies and explained the pros and cons of each 
one, supported by microwave simulation techniques that are specifically suited to deal 
with dimensions shrinking from millimetre to nanometre scale. We have then described 
a new structure that combines a good control of microwave transmission modes with a 
maximized value of the magnetic field available to drive coherent control of a spin qubit. 
This structure has been employed to achieve coherent control of a single-donor electron 
spin qubit in silicon. The experiments show that the simulation describes accurately 
the magnetic field produced by the microwave loop, whereas the electric field modelling 
remains in need of better understanding. 

We expect the design guidelines presented in this paper to be useful in assisting 
and facilitating the intense efforts towards performing coherent spin control in 
nanostructures, where high values of oscillating magnetic fields at frequency 3> 10 GHz 
must coexist with charge-sensitive devices. 
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